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ABSTRACT 

The disruption of a binary star system by the massive black hole at the Galactic 
Centre, SgrA*, can lead to the capture of one star around SgrA* and the ejection of 
its companion as a hypervelocity star (HVS). We consider the possibility that these 
stars may have planets and study the dynamics of these planets. Using a direct A-body 
integration code, we simulated a large number of different binary orbits around SgrA*. 
For some orbital parameters, a planet is ejected at a high speed. In other instances, 
a HVS is ejected with one or more planets orbiting around it. In these cases, it may 
be possible to observe the planet as it transits the face of the star. A planet may also 
collide with its host star. In such cases the atmosphere of the star will be enriched 
with metals. In other cases, a planet is tidally disrupted by SgrA*, leading to a bright 
flare. 
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1 INTRODUCTION 



Hypervelocit y stars (HVSs) were first theorized in 1988 
Hills 1 19881). and discovered observationally in 2005 



Brown et al. i l2005|V At least 16 HVSs have been identified 



Brown et al. 


2006a; iBrown et al. 2006bl; iBrown et al.l 


2007; 


Brown et al. 


2009). There are a number of proposed mech- 



anisms for the production of HVSs. The best studied and 
arguably most likely mechanism is the Hills' scenario where 
the close interaction of a binary star system and massive 
black hole (MBH) can produce a HVS with sufficient ve- 
locity to es c ape the gravitational p ull of the Milky Way 
jHillsl Il988l; iGould fc Quilled 120031; [Yu fe Tremaine 12003 ; 
Ginsburg fc Loebl |200d; IPerets et al) 120071: iPerets! (2009a 



Madigan et al.ll2009l . l201ll ). Other mechanisms include the 



inter action of stars with stellar black holes |0'Learv fc Loebl 
120081). the inspiral of an in t ermediate mass black hole 
jHansen fc Milosavljevklj2003l ; lYu fc Tremaine! 120031 ; iLevinl 

and the disrupt i on of a triple 



Sesana et al 



I2OO6I ; ... 

systeniTLu et al.l 2007|: 



I Ginsburg fc Peret 



2009) 



Sesana et al. 2009; Pcrcts 2009b; 



20111 "). We focus on the disruption of a 



tightly bound binar y by a central black hole of inferred mass 
~4x 10 6 M (e.g. ISchodel et al]|2003l ; iReid fc Brunthalerl 
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| 2004[ ; iGhez et all 120051 ; iGhez et al l 120081 ; iGillessen et al] 
120091 ). 

Simulations show that tight binaries disrupted by the 
MBH can result in the ejection of one component with ve- 
locities comparable to those of the observ ed HVSs (e.g. 
iGinsburg fc Loebl 120061 ; iBromlev etai]|2006D. and perhaps 
significantly larger velocities as well ( Sari et al]|2010l ). When 
a hypervelocity star is produced, the companion to the 
HVS is left in a high ly eccentric orbit around SgrA* 
(|Ginsburg fc Loebl |2006 ) . Furthermore, a small fraction of 
binary systems may collide when disrupted, and if the 
collision velocity i s sma ll enough the system can coalesce 
|Ginsburg fc Loebl 120071 : lAntonini et all l201ll ). For plan- 
ets orbiting a star, tidal dissipation will eventually result 
in t he infall of short- p eriod planets into t heir host star 
(e.g. iRasio et~al 1 1 19961 ; iJackson et all 120091 ); however the 
timescale is on the order of Gyrs, and planets with orbital 
periods of < 1 day are predic ted to survive while the host 
star is on the main-sequence (|Hansenll2O10l ). A number of 
Jupiter-mass planets have already been found with orbita l 
periods < 1 day (e.g. iHebb et~alj|2010l ; iHellier et al.ll201ll ). 
Thus, it is only natural to consider the scenario where plan- 
ets are orbiting such tightly bound binary systems which are 
subsequently disrupted by the MBH. Our goal is to examine 
the possible orbital dynamics and determine whether a HVS 
may host a planetary system. 

In §2 we describe the codes and simulation parameters. 
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In §3 we discuss the origin of hypervelocity planets (HVPs) 
and transits around HVSs, and in §4 we discuss some pos- 
sible outcomes for the disruption of a binary system with 
planets. Our goal was not to cover the entire available phase 
space, but to determine whether some tight binaries with 
planets could produce transits or other observable effects. 



2 COMPUTATIONAL METHOD 

In our study we have used the publically available N-body 
code written bv lAarsetbl l|l999l ) n. We adopted a small value 
of 1(T 8 for the accuracy parameter 77, which determines the 
integration step through the relation dt — \J r\F / (d 2 F / dt 2 ) 
where dt is the timestep and F is the force. The soften- 
ing parameter, eps2, which is used to create the softened 
point-mass potential, was set to zero. We treat the stars 
and planets as point particles and ignore tidal and general 
relativistic effects on their orbits, since these effects are small 
at the distance (> 10AU) where the binary star system is 
tidally disrupted by the MBH. In total, we ran over 8000 
simulations. 

We set the mass of the MBH to M = 4x 1O 6 M , and the 
mass of each star to m = 3Mf), which is comparable to the 
known masses of HVSs (|Fuentes et al.ll2006l ; iPrzvbilla et ail 
2008). Since the planet's mass <C m* we treat all planets as 
test particles. All runs start with the centre of the binary 
system located 2000 AU (= 10~ 2 pc) away from the MBH 
along the positive y-axis. This distance is larger than the 
binary size or the distance of closest approach necessary to 
obtain the relevant ejection velocity of HVSs, making the 
simulated orbits nearly parabolic. Our simulations include 
orbits in a single plane as well as simulations with the binary 
coming out of the orbital plane at 90°. For simplicity, the 
planetary orbital plane was kept the same as that of the 
binary. We used the same initial distance for all runs to make 
the comparison easier to interpret as we varied the distance 
of closest approach to the MBH or the relative positions of 
the stars and planets within the system. 

We ran two primary sets of simulations. The first set 
consisted of a binary system with two planets. Each star had 
one planet initially at a p = 0.02 AU from its host star on a 
circular orbit (with eccentricity e = 0). The second set con- 
sisted of a binary system with four planets. Initial conditions 
were similar to the previous data set, with the exception 
that the second planet was placed on a circular orbit with 
distance a p = 0.03 AU from the star. We varied the initial 
separation between the two stars from a* = 0.05 to 0.5 AU 
and precluded tighter binaries for which two stars develop a 
common envelope and coale sce. Similarly, much wider bina- 
ries may not produce H VSs (|Hillslll99ll ; lBromlev et al.ll2006l : 
iGinsburg fe Loebl l2006) . However, in our simulations we as- 
sume planetary distances significantly shorter than 0.02 AU 
would bring the planet into the star's atmosphere and thus 
are precluded. Similarly, too large a distance could result in 
a dynamical instability whereby the planets are ejected from 
the system or collide with a star. In order to determine the 
effect of the planetary separation on the orbital dynamics, 
we ran additional simulations with a fixed binary semimajor 



1 http: / /www. ast.cam.ac.uk/~sverre/web/pages/nbody.htm 



axis a* = 0.2 AU and varied the planetary separation in the 
range a v = 0.02-0.06 AU. Our simulations excluded planets 
around single stars since such systems will never produce 
HVSs. 

A three-body encounter can lead to chaos, and the ini- 
tial phase of the stellar orbit can greatly vary the outcome 
I Ginsburg fc Lo cb 2006). Therefore, we sampled cases with 
initial phase values of 0-360° at increments of 15°. We gave 
the binary system no radial velocity but a tangential velocity 
with an amplitude such that the effective impact parameter 
(b) is in the range of 5-35 AU. The distribution of impact 
parameters is as follows: our simulations had 6 = 5, 10, 15, 20 
AU with a likelihood of ~ 20% of all runs for each value, 
while the remainder at 6 = 25 or 30 AU had a ~ 10% likeli- 
hood each. We expect no HVSs to be produced at subst an- 
tially larger impact parameters l|Ginsburg fc Loebll2006l ). 



3 ORIGIN OF HYPERVELOCITY STARS AND 
HYPERVELOCITY PLANETS 

Given a binary system with stars of equal mass m separated 
by a distance a and a massive black hole (MBH) of mass 
M S> m at a distance 6 from the binary, tidal disruption 
would occur if b < bt where 



M 



(1) 



The distance of closest approach in the initial plunge of the 
binary towards the MBH can be obtained by angular mo- 
mentum conservation from its initial transverse speed v± at 
its initial distance from the MBH, d, 



v±d ■■ 



1/2 



(2) 



The binary will be tidally disrupted if its initial transverse 
speed is lower than some critical value, 



v± 



< 



(GMa 



,1/2 



d 



M 



1/6 



-1/2 



= 10 



1/6 T 
V 5 "3.3 



km s 



(3) 

where a_i = (a/0.1 AU), d 3 . 3 = (d/2000 AU), m . 5 = 
(m/3M Q ), and we have adopted M = 4 x 10 6 M Q . For 
M -L U -L,crit, one of the stars receives sufficient kinetic en- 
ergy to become unbound, while the second star is kicked into 
a tighter orbit around the MBH. The ejection speed, v c j, of 
the unbound star can be obtained by considering the change 
in its kinetic energy ~ vSv as it acquires a velocity shift of 
order the binary orbital speed Sv ~ \J Gm/a during the 
disruption process of the binary at a distance ~ b t from the 
MBH when the binary centre-of-m ass speed is v ~ v/ GM/bt 
<|Hillslll98Sl ; lYu fc Tremainell2003l '). At later times, the binary 
stars separate and move independently relative to the MBH, 
each with its own orbital energy. For v < «j_, cr it, we therefore 
expect 



(?r(f) 



1/2 



1/2 



1.7 x 10 3 m^V 



1/2 
1 



km s 



(4) 

A planet with mass m p <C m* may also be ejected at 
even higher speeds. The mechanism that ejects HVPs in- 
volves the interaction of the MBH, the binary star system, 
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Velocity Distribution 




v (km s 



Figure 1. Velocity distribution of HVSs and HVPs. This sample 
comes from 1000 simulations. The initial binary separation for 
each system was a* =0.2 AU, and we let the planetary separation 
vary with uniform probability in the range a p = 0.02-0.04 AU. 
These runs are strictly for systems with two planets, however we 
get similar distributions for systems with four planets. The lowest 
HVP velocities are ~ 700 km s — 1 , which corresponds with the 
lowest velocities for HVSs. The average HVP velocity is ~ 3000 
km s -1 , and the average HVS velocity is ~ 1500 km s (see also 
Table [TJ. The HVPs are denoted in light blue, and the HVSs in 
the darker color. The overall shape of the distribution is similar in 
both cases, and in both cases there are outliers. Note that there 
are ~ 3 HVPs for each HVS. 



and the planet, and does not admit a simple analytical es- 
timate. Based on our simulations, the velocities of hyper- 
velocity planets (HVPs) are on average ~ 1.5-4 times the 
velocity of a typical HVS (see Table [IJ. We found a few ex- 
amples of HVPs with exceptionally high speeds, v ~ 10 4 km 
s _1 , however as observed in Figure [1] these are rare with a 
probability < 1%. Although Figure [T] represents one specific 
parameter set, we do get comparable results for other runs. 



4 THE FATE OF PLANETARY SYSTEMS 

We analyze statistically the orbital properties of the binary 
star system after it is disrupted by the MBH. We ran sim- 
ulations both planar and with the binary out of the orbital 
plane at 90°. In our simulations the initial planetary semi- 
major axis is constrained. Too small a separation would lead 
to a plunge of the planet into the host star's atmosphere, and 
too large a separation would result in an unstable configu- 
ration in which the planet ultimately collides with a star or 
else is ejected from the system. 

Figure [5] illustrates two possible outcomes. In both in- 
stances, HVSs with orbiting planets are produced (dashed 
line), while the companion star (solid line) orbits the MBH 
in a highly eccentric orbit. The companion star's planets are 
themselves removed, and either stay bound to the MBH in 
different orbits or are ejected as HVPs. Table [2] shows vari- 
ous outcomes from our simulations including the fraction of 
HVPs, planets around HVSs, and free planets around the 
MBH. 



a* (AU) 


VHVS (km s 1 ) 


"U7P (km s 1 ) 


0.05 


2700 


3800 


0.10 


2000 


3100 


0.20 


1400 


3500 


0.30 


1300 


4100 


0.40 


1400 


4100 


0.50 


1100 


4400 


0.05 






0.10 






0.20 


1500 


3300 


0.30 


1300 


3900 


0.40 


1400 


4200 


0.50 


1100 


4200 



Table 1. Average velocity of HVSs (second column) and HVPs 
(third column) for different values of a*. The top four rows show 
the values obtained from our simulations with two planets, and 
the bottom rows our results with four planets. For our simulations 
with four planets, the outer planets either collided with a star or 
were immediately ejected when a* < 0.2 AU, hence no values are 
listed for 0.05 and 0.10 AU. All values are uncertain to within 
±50 km s" 1 



4.1 HVPs 

Our simulations show that for tight binaries with semimajor 
axes in the range a* = 0.05-0.5 AU and planets with semi- 
major axes a p = 0.2-0.05 AU, the probability of producing 
HVPs is high (see Table |2J. For a binary with two planets, 
we find that the probability of producing a HVP is ~ 30- 
40%. For a binary with four planets, the probability rises to 
~ 70-80%. However, since planets are extremely faint, it is 
impossible to detect HVPs directly with existing telescopes. 



4.2 Transits Around HVSs 

The probability of observing an eclipse as a planet orbits its 
host star is the transit probability (Pt) given by 



Pt = 



Ri, ± Rp\ /l + esinw 
1 -e 2 



(5) 



where J?* is the star's radius, Rp the planet's radius, e 
is the eccentricity, and u) is the arg ument of periapse (see 
iMurrav fc Correial[201ol : IWinnll201of l. Note that the + sign 
in equation ((5]) includes grazing eclipses, which are otherwise 
excluded. When R+ 2> Rp and e = 0, we get 



0.015 



/ R* 



1AU 



(6) 



For massive (m IMq) stars on the main- sequence, it 
is safe to assume R* S> Rp and thus the probability of 
a transit is independent of planetary radius. However, the 
fractional flux decrement is given by (Rp/R*) 2 . Therefore a 
larger planetary radius is important in order to be able to get 
proper photometry. We assume a Jupiter-sized planet, thus 
(Rp/R+) 2 ~ 1%. Our simulations produced planets orbiting 
HVSs with e = 0.6. Although an eccentricity will increase the 
probability of a transit, equation (0 is sensitive to high ec- 
centricity and for a limiting case it is reasonable to use e = 0. 
The timescale for eccentricity damping, whereby a planet's 



© 0000 RAS, MNRAS 000, 000-000 



4 Ginsburg, Loeb, Wegner 




Figure 2. The two panels illustrate possible outcomes after a binary system with planets is disrupted by the MBH. The MBH is located 
at the origin, and the binary system starts at an initial distance of 2000 AU along the positive y-axis. In both cases the initial binary 
separation was a* = 0.2 AU and the planetary separation a p = 0.02 AU for the innermost planets and a p = 0.03 AU for the outer planets. 
Left: After binary disruption a HVS is produced with two bound planets, as marked by the dashed line. The second star, marked by 
the solid line, stays in a highly eccentric orbit around the MBH. The second star's planets are removed, and the first falls into a highly 
eccentric orbit close to the MBH, while the second is ejected into a much larger, but also highly eccentric orbit around the MBH. Right: 
After binary disruption a HVS is produced with two planets in orbit as marked by the dashed line. The second star, marked by the solid 
line, remains in a highly eccentric orbit around the MBH. The second star's planets are ejected as HVPs. 



a* (AU) 


HVSs 


HVPs 


HVSs + Planets 


Bound Stars + Planets 


Free Planets 


0.05 


0.40i0.03 


0.36±0.03 


0.05±0.01 


0.08i0.01 


0.80i0.03 


0.10 


0.37±0.03 


0.36±0.03 


0.08±0.01 


0.35i0.03 


0.63i0.03 


0.20 


0.20±0.02 


0.42±0.03 


0.05±0.01 


0.40i0.03 


0.58i0.02 


0.30 


0.07±0.01 


0.39±0.03 


0.02±0.01 


0.35i0.03 


0.62i0.03 


0.40 


0.03±0.01 


0.35±0.03 


0.01i0.01 


0.39i0.03 


0.62i0.03 


0.50 


0.02±0.01 


0.43±0.03 


0.01i0.01 


0.42i0.03 


0.57i0.02 


0.05 












0.10 












0.20 


0.19i0.02 


0.83±0.04 


0.11i0.02 


0.59i0.04 


0.64i0.02 


0.30 


0.05±0.01 


0.70±0.04 


0.03±0.01 


0.59i0.04 


0.67i0.02 


0.40 


0.02±0.01 


0.71±0.04 


O.OliO.Ol 


0.63i0.04 


0.65i0.02 


0.50 


0.02±0.01 


0.78±0.04 


O.OliO.Ol 


0.76i0.04 


0.60i0.02 



Table 2. Probability for various outcomes in our ensemble of runs. The top rows show the results for simulations with two planets, and 
the bottom rows for our simulations with four planets. Note that for our simulations with four planets, the outer planets either collided 
with a star or were immediately ejected when a* < 0.2 AU, hence no values are listed for 0.05 and 0.10 AU. The distance between each 
planet and host star is a v = 0.02 AU, and for simulations with four planets the second set of planets had a p = 0.03 AU. The first column 
is the initial distance between the two stars. The second column shows the fraction of HVSs produced, and as expected tighter binaries 
produce more HVSs. The third column shows the rate of HVP production. The fourth column shows the probability of producing HVSs 
with planets in orbit. This column is also represented graphically in the middle panel of Figure [3] The fifth column shows the fraction 
of stars that are bound to the MBH and have planets in orbit which may produce transits (see Section l4.4l l. The last column shows the 
fraction of free planets around the MBH. Probabilities in each row are independent of each other. Note that the sum of the probabilities 
for each row add to more than one due to the fact that we have two or four stars for each simulation. Therefore, as illustrated in Figure 
[2] there may be multiple outcomes for any given simulation. The quoted values include Poisson errors. 



eccentricity is circularized is not well understood, and de- 
pending on initial conditions one can get either short (~ 10 
Myr) or long (~ Gyr) circularization times l|Naozl [20121 ). 
Furthermore, e = is consistent with the observed eccen- 
tricities for most extrasolar planets with a v < 0.1 AU (e. g. 
Ijackson et al.|[200 1; lHansenll2oTol ; iMatsumura et al.ll2010h - 

Table [2] shows the probability of producing a HVS with 



orbiting planets. This probability is calculated as the frac- 
tion of all the runs we conducted in the restricted range of 
parameters for which HVSs are produced. We used a bi- 
nary semimajor axis between a* = 0.05-0.5 AU. Smaller a* 
would lead to a common envelope and are precluded. Larger 
a* are likely to either eject the planet from the system be- 
fore producing a HVS, lead to a collision, or fail to produce 
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a HVS altogether. The results for our simulations with two 
planets and four planets are qualitatively similar with some 
expected statistical variations, as illustrated in the middle 
panel of Figure [3] Our simulations also indicate that the 
planetary semimajor axis needs to be between a p = 0.02- 
0.05 AU in order for a planet to remain bound to a HVS. 
Smaller a p are precluded due to the fact that they would sink 
into the star's atmosphere, and at larger a p the planets are 
tidally removed before the HVS is produced. For two-planet 
systems with a* = 0.2 AU, we find that when a p = 0.02- 
0.05 AU the probability of producing a HVS with planets is 
constant. We get similar results for four planets, however the 
probability starts to rapidly decrease when a p = 0.04 AU. 
These results are illustrated in the right panel of Figure [3] 
In total, we find that the probability of producing a HVS 
with planets peaks at ~ 10%. Assuming a HVS has planets 
in orbit, the probability of observing a transit is given by 
equation We plot these probabilities in Figure 2] 



star's surface. The metallicity of stars with planetary sys- 
tems has been shown to be greater than expect ed by a fac- 
tor of ~ 2 fe.g. lSantos et al.ll2003l ; [Vauclairll2004h . Such high 
metallicity may be due to tidal interacti ons that have caused 
exoplanets to fall in to their host stars jLevrard et al.l [20091 : 
I Jackson et al.| [2009). Although the mechanism by which a 
star may lose its high- m e tallicity materi al is not well un- 
derstood (|Sasselovll20l2T ), iGaraudl (120111 ) has shown that 
the relative metallicity enhancement on the surface due to 
the infall of exoplanets, survives a much shorter period for 
higher-mass stars. A star of 1.5M© would lose ~ 90% of its 
enhanced metallicity in merely 6 Myr. Thus, a HVS (typi- 
cally found at large distances from their origin at the Galac- 
tic Centre) may have lost nearly all its enhanced metallic- 
ity relatively ea rly. Furthermore, meta l-rich HVSs have not 
been detected (IKo llmeier et al.l l201lh ■ We suggest that if 
current or future HVSs are found with enhanced metallic- 
ity, it may be due to a relatively recent planetary collision. 



4.3 Collisions 



4.4 Bound and Free Planets Around SgrA* 



Assuming that a binary is disrupted at a random angle and 
ignoring gravitational focusing, the probability for a collision 



4R* 



27Ti? s , 



(7) 



where i?* is the radius of the star (assuming radii are ap- 
proximately equal) and R sep is the average separation of the 
stars. The stars will merge if the relative velocity of impact 
is less than the escape velocity from the surface of the star 
(~ 500 km s _1 ). A simple estimation for the impact velocity 
comes from conservation of energy, 



1 m 1 m 2 .2 Gmim 2 

h — r = const, 

2 mi + mi r 

which yields the relative impact velocity for two stars 



v f = 



2G(mi +m 2 ) 



1 



1 



1/2 



Similarly, for a planet around a star, we arrive at 



Vf 



Gm t 



2a 



1 1/2 



(8) 



(9) 



(10) 



Equations © and (|10p yield a velocity similar to the escape 
speed from the surface of the star v ~ 500 km s" 1 , however 
they do not take tidal forces nor gravitational focusing into 
account, therefore the actual impact velocities are expected 
to be slightly higher. For a Jupiter-like planet, a direct col- 
lision will release at minimun 10 45 erg, and thus might be 
observable as a flare. 

The rate of HVS productio n in the Milky Way is 
~ 10~ 5 yr" 1 (|Brown et alj|2006bl). and the total co llisional 
rate between stars is ~ 10 % (|Ginsburg fc Loebll2007l ). Based 
on our simulations, we find that the collisional probability 
for a system with two planets is ~ 0.1%, whereas for a sys- 
tem with four planets it is ~ 1%. For systems with four 
planets, the collisions are almost entirely due to the two 
outer planets. 

Should a planet collide with its host star, the planet 
will be destroyed and add high-metallicity material to the 



When a binary system is tidally disrupted by the MBH, a 
HVS may be produced while the companion star remains in 
a highly eccentric orbit around the MBH. If the disruption 
fails to produce a HVS, both stars will orbit the MBH. In 
either case, our simulations show a substantial probability 
for a planet to remain bound to a star that is in orbit around 
SgrA* (see Table • It is unclear whether in the long run 
such planets will be dynamically stable or eventually fall into 
their host stars. However, it is certainly possible that some of 
the stars near SgrA* host planetary systems. Should a star 
near the Galactic Center have a planet transit the surface, 
the change in magnitude will be 8m ~ 0.01. Near-infrared 
photometry can get close to 0.01 mag precision, however 
extinction an d the crowding of st ars make such observations 
difficult fe.g. lSchodel et alj|2010l ). 

The probability for at least one planet to fall into a 
highly eccentric orbit around SgrA* is on average > 60% for 
all our runs. If such a planet passes within the tidal radius 
of the MBH 



Rt ~ Rr> 



M A 



1/3 



(11) 



(where R p and m p are respectively the radius and mass of 
the planet), the tidal force will disrupt the planet. A star 
will similarly be disrupted if it passes too close to its tidal 
radius, and might produce an optical flare of ~ 10 43 erg s _1 
jStrubbe fc Quataerdlioogl ). For a planet disrupted by the 
MBH, the accretion rate M mass would be lower, leading to 
a luminosity 

L = eMc, (12) 

where e is the radiative efficiency. IZubovas et al.l (|201 ll ) sug- 
gest that asteroids disrupted by SgrA* may be responsible 
for flares with L ~ 10 34 - 10 39 erg s" 1 . 
note that a brighter flare (L ~ 10 41 erg s 
to have occurred ~ 300 years ago may be due to the tidal 
disruption of a planet. Figure [5] shows the probability for 
a planet in orbit around the MBH to be tidally disrupted. 
The probability for a system with two planets is shown by 
the dashed line, and four planets by the solid line. 



Furthermore, they 
_1 ) that is inferred 
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Probability of a HVS with Planets vs the Impact Parameter Probability of a HVS with Planets vs Binary Separation Probability of a HVS with Planets vs Planetary Separation 




Figure 3. Probability of producing a HVS with orbiting planets versus various parameters. Left: Probability as a function of impact 
parameter, b, relative to SgrA*. This cumulative plot shows that an impact parameter between 10 — 20 AU is optimal. The probabilities 
here are defined relative to all our HVS runs. Middle: Probability as a function of the initial distance between the two stars. In this 
scenario we exclude simulations with four planets when a* < 0.2 AU, due to the fact that at such distances the two outermost planets 
were disrupted and either collided with a star or were ejected. The initial planetary separation is a p = 0.02 AU for all runs with two 
planets, and for our simulations with four planets the second set of planets had initial distance a p = 0.03 AU. Right: Probability as 
a function of initial planetary separation, a p . We precluded orbits with a p < 0.02 AU, since such orbits would take the planet into the 
star's atmosphere. The results for each figure are averaged for all orbits, those in a single plane and orbits out of the plane. The values 
of a p displayed along the x-axis correspond to the first set of planets. For simulations with four planets, the second set was placed 0.01 
AU from the first. Note that for the middle and right panels the impact parameter is not constrained. 



Probability of a Transit vs Binary Separation Probability of a Transit vs Planetary Semimajor Axis 




a.(AU) a p (AU) 



Figure 4. Probability of observing a transit for a planet orbiting a HVS. Left: Probability of a transit (log Pt) vs the initial binary 
separation (a*). We exclude simulations with four planets when a* < 0.2 AU, since the outermost planets are unstable at these parameters. 
Right: Probability of a transit (log Py) vs the initial distance between the planet and host star (a p ). In this scenario a* is fixed at 0.2 
AU. In both instances, our results for two planets and four planets are qualitatively similar. The probability of observing a transit ranges 
from 0.35 - 0.90. 



5 CONCLUSIONS 

Our simulations show that in order to produce a HVS with 
planets in orbit, the initial binary separation needs to be in 
the range a* = 0.05-0.5 AU, and the planetary separation in 
the range a p = 0.02-0.05 AU (e.g. "hot Jupiters"). For such 
parameters there is up to a ~ 10% probability that a HVS is 
produced with orbiting planets (see Table[5J). In such scenar- 
ios there is a high probability that at least one planet tran- 
sits the star (see Figure [4]). For hot Jupiters, the fractional 
flux decrement in the light curve, (R P /R*) 2 ~ 1%, would 
typically last for hours. Assuming no systematic errors, a 
powerful enough telescope can get milli-magnitude photom- 



etry (corresponding to fractional flux sensitivity of 0.1%) 
and hence determine whether or not a HVS has any tran- 
siting planets. Our examples of planets around HVSs have 
orbital periods in the range 0.3 to 19 days, requiring observ- 
ing programs of this duration. The Hills' mechanism that 
produces such HVSs also produces HVPs. Unfortunately, it 
is extremely difficult to detect a free HVP. Our simulations 
indicate that HVPs may achieve speeds ~ 10 4 km s _1 in 
rare circumstances. 

When a binary system with our given parameters is 
disrupted, it is possible that a planet will collide with its 
host star. In such instances the surface may be enriched with 
metals. Such an enrichment is not believed to last more than 
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Figure 5. Probability for a planet freed during the production of a HVS from a binary star system, to be tidally disrupted by the MBH. 
Left: Probability as a function of the binary separation. We exclude simulations with four planets when a* < 0.2 AU due to instability. 
Right: Probability as a function of planetary separation. a p = 0.02 AU was the minimum semimajor axis length. For a v = 0.06 AU 
the planets either crashed into a star or were removed, hence the probability is ~ 0. The origin of the break between the simulations 
with two planets (dashed line) and four planets (solid line) is due primarily to the fact that for four planets, the outer planets are easily 
disrupted with increasing a*, and are ejected or crash into a star before being tidally disrupted by the MBH. For increasing a p , both 
two planet and four planet simulations have similar stability, however there are more planets in the latter case and thus the probability 
does not drop as quickly. 



a few Myr for massive stars |Garaudll201ll ). The age of HVSs 
exceed the time necessary for the excess metallicity in the 
star's atmosphere to be erased. If observations show that 
a HVS has an unusually metal-rich atmosphere, this may 
indicate that a planet has relatively recently fallen into the 
star's atmosphere. 

Planets that are not ejected as HVPs and are not or- 
biting any HVSs, may still be detectable as transits around 
stars orbiting SgrA*. Assuming a Jupiter-like planet, the 
change in magnitude will be 5m ~ 0.01. There is also 
a high probability that free-floating planets will be pro- 
duced. It may not be possible to ascertain whether a free- 
floating planet was originally bound to a binary system 
that was disrupted. However, our simulations show that 
free-floating planets from disrupted systems should be in 
highly eccentric orbits around the MBH with an eccentric- 
ity ~ 0.96 or greater, sim ilar to former companions to HVSs 
jGinsburg fc Loebll2006h . Furthermore, it is also possible for 
a planet or star to be tidally disrupted by the MBH, pro- 
ducing a substantial flare from SgrA*. The luminosity and 
duration of the flare would indicate whether it was a planet 
(L ~ 10 40 -10 41 erg s" 1 ) or star (L ~ 10 43 erg s" 1 ) that fe d 
the MBH (|Strubbe fc Quataertj[2009l ; IZubovas et al.ll201ll ). 

The detection of even one planet around a HVS can 
shed light on planetary formation and evolution within the 
central arcsecond of SgrA*. In particular, it would support 
the notion that young stars near the MBH have close-in 
planets. It has been suggested that planets and even aster- 
oids may form in protoplanetar y discs around stars orbiting 
a MBH (|Navakshin et alll2012h . Rec ently, a gas cloud wa s 
discovered plunging towards SgrA* (|Gillessen et al.ll2012t ). 
While it was suggested that this clo ud originated from s tellar 
winds or from a planetary nebula (jBurkert et al.ll2012r ). an- 
other theory is that it originated from a proto-planetary disk 



around a low-mass star, which implies that planets may form 
in the Galactic Centre (jMurrav-Clav fc Loebll2012r i. The de- 
tection of a planet transiting a HVS would lend considerable 
support to this interpretation. If most stars in the Galactic 
centre form with planetary systems, then our simulations in- 
dicate that there should be planets in highly eccentric orbits 
around SgrA*. 
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